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We investigate cosmological behavior in the quasi-Dilaton non-linear massive gravity. We perform 
a detailed dynamical analysis and examine the stable late-time solutions relevant to late time cosmic 
acceleration. We demonstrate that a dark-energy dominated, cosmological-constant-like solution is 
a late attractor of the dynamics. We also analyze the evolution of the universe at intermediate times, 
showing that the observed epoch sequence can be easily obtained in the model under consideration. 
Furthermore, we study the non-singular bounce and the cosmological turnaround which can be 
realized for a region of the parameter space. Last but not least, we use observational data from Type 
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I. INTRODUCTION 

It is widely believed that an elementary particle 
of mass 77i and spin s is described by a field which 
transforms according to a particular representation of 
Poincare group. In field theory, formulated in Minkowski 
space-time, mass can either be introduced by hand or 
via spontaneous symmetry breaking. The spin-2 field 
hfu, is supposed to be relevant to gravity as it shares 
an important property of universality with Einstein gen- 
eral relativity (GR) a la Weinberg theorem and GR can 
be thought as an interacting theory of h^ v . It is there- 
fore natural to first formulate the field theory of h^ w in 
Minkowski space-time and then extend the concept to 
non-linear background. 

The first formulation of a linear massive theory of grav- 
ity was given by Fierz and Pauli |l| with a motivation to 
write down the consistent relativistic equations for higher 
spin fields including spin-2 field. The linear theory of 
Fierz-Pauli, however, suffers from van Dam, Veltman, 
Zakharov (vDVZ) discontinuity @,y|: in the zero-mass 
limit, the theory is at finite difference from GR. The un- 
derlying reason of the (vDVZ) discontinuity is related to 
the fact the longitudinal degree of freedom of massive 
graviton does not decouple in the said limit; it is rather 
coupled to the source with the strength of the universal 
coupling at par with the massless mode. 

It was pointed out by Vainshtein that linear approxi- 
mation in the neighborhood of a massive source breaks 
down below certain distance dubbed Vainshtein radius 
and by incorporating non-linearities one could remove the 
discontinuity [J] present in the linear theory. However, 
the sixth degree of freedom, suppressed in Pauli-Fierz 
theory, which is essentially a ghost known as Boulwarc- 
Deser (BD) ghost [fj, becomes alive in the non- linear 



background. 

Efforts were then made to find out the non-linear ana- 
log of Fierz-Pauli mass term requiring the ghost to be 
systematically removed. The authors of Ref.@, 0] dis- 
covered a specific nonlinear extension of massive gravity 
dRGT) that is BD-ghost free in the decoupling limit(scc 
8] for a review). Subsequently, it was demonstrated 
9( that the Hamiltonian constraint is maintained at the 
non-linear order along with the associated secondary con- 
straint, which implies the absence of the BD ghost. Apart 
from the field theoretic interests of such a construction, 
there is a cosmological motivation to study massive grav- 
ity. This formulation provides a new class of (Infra-Red) 
gravity modification which can give rise to late-time cos- 
mic acceleration and it is therefore not surprising that 
massive gravity has recently generated enormous inter- 
est in cosmology |9l-r71j. 

In spite of the the successes of the framework, it 
has challenging problems to be resolved, namely, the 
model does not admit spatially fiat FLRW background 
whereas isotropic and homogeneous cosmolo gica l solu- 
tions (K = ±1) are perturbatively unstable [72j |. The 
latter speaks of some inherent difficulty of dRGT. This 
led the authors of [73, uM t° abandon isotropy whereas in 
Ref . [75| , effort was made to address the problem by mak- 
ing the graviton mass a field dependent quantity through 
a more radical approach, namely by extending the the- 
ory to a varying graviton mass driven by a scalar field 
[7a . |77| . A similar approach was followed in [78[ , where 
the framework of massive gravity was extended by in- 
troducing a quasi-dilaton field. It was pointed that the 
model could rise to a healthy cosmology in the FLRW 



background 1 . 

In this work we carry out detailed dynamical investi- 
gations of cosmological behaviors in the aforementioned 
quasi-Dilaton non-linear massive gravity. In particular, 
we perform dynamical analysis relevant to late time cos- 
mic acceleration, and moreover we examine the realiza- 
tion of non-singular bouncing solutions. Finally, we use 
observational data in order to constrain the model pa- 
rameters. 

The plan of the paper is as follows: In section [TT| we 
review the quasi-dilaton non-linear massive gravity and 
its cosmological implications. In section ITTTl we perform 
a detailed dynamical analysis and obtain the late-time 
stable solutions, while in section IIVI we investigate the 
universe evolution focusing on bouncing and turnaround 
solutions. In section [V] we use observations in order to 
constrain the model parameters and finally in section |VT1 
we summarize the results of our analysis. 



the action corresponding to non-linear massive term, 
along with the action of a massless dilaton a 



S a = -^[d*x(daf 



(4) 



In action (|3|), as, an are two arbitrary parameters and 
Ui are specific polynomials of the matrix 



K£ = 5» - e a ' Mp > V9^ a d a rd^ b Vab, (5) 

given by 

U 2 = 4([/C] 2 - [K 2 ]) (6) 

[/ 3 = [/C] 3 -3[/C][/C 2 ]+2[/C 3 ] (7) 

U A = [/C] 4 - 6[/C] 2 [/C 2 ] + 3[/C 2 ] 2 + 8[/C][/C 3 ] - 6[/C 4 ]. (8) 

In ©, rial, (Minkowski metric) is a fiducial metric, and 
4> a (x) are the Stiickelberg scalars introduced to restore 
general covariance [8l[ . 



II. QUASI-DILATON NON-LINEAR MASSIVE 
GRAVITY AND COSMOLOGY 



B. Cosmology 



In this section we briefly review the quasi-dilaton non- 
linear massive gravitational formulation following [78| . 
In the first subsection we present the gravitational the- 
ory itself, while in the next subsection we examine its 
cosmological implications in presence of matter and ra- 
diation. 



A. Massive gravity with quasi-dilaton 

Let us consider the following action of non-linear mas- 
sive gravity with dilaton 



iSer — Seh + <S n 



<-><Tj 



with 



'EH 
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I d^xy/^g R 



the usual Einstcin-Hilbcrt action, 
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(i) 



(2) 



(3) 



1 During the finalization of the present work, a perturbation anal- 
ysis of the stability of the model was performed in [79l. ISOtj . where 
it was shown that the de Sitter solution of the model may have 
a ghost instability for short wavelength modes, shorter than cos- 
mological scales. The calculations were performed at the linear 
order of perturbations, or equivalcntly at second order in the ac- 
tion. It will be interesting to see whether the ghost exists in a 
non-linear approach of the stability analysis. 



In order to apply the above gravitational theory in cos- 
mological frameworks, we also have to incorporate the 
matter and radiations sectors. Thus, the total action we 
consider is as follows 



S = S } 



EH 



Sn 



O ct + S m + S r , 



(9) 



in which S m and S r denote the standard matter and ra- 
diation actions, corresponding to ideal fluids with energy 
densities p m , p r and pressures p m , p r respectively. In 
summary, we consider a scenario of usual nonlinear mas- 
sive gravity [f| 0] coupled with a dilaton field a, where 
the coupling is introduced through (|5|) [78 1. 

Next, we consider a flat Friedmann-Lemaitre- 
Robcrtson- Walker (FLRW) metric of the form 

ds 2 = g^dx^dx" = -N{tfdt 2 + a 2 (t)6 ij dx i dx :i ) (10) 

while for the Stiickelberg scalars we consider the ansatz 

0° = /(*), f = x\ (11) 

In this case the actions ©-© become: 

Seh = -3Mp, f dt 



aa 



*->mass — 3m Mpi 
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dt a 3 
where we have defined 
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(12) 
(13) 
(14) 



Gi(0 = (1 - 2 - £ 4 ^(1 - 0(4 - + "4(1 - 2 



l + i a 3 (l-0 + «4(l-0 2 



(15) 
(16) 



being functions of the introduced dilaton- variable 



£ 



o<r/M F 



(17) 



Variation with respect to / gives the constraint equa- 
tion (as an essential feature of any model of non-linear 
massive gravity) 



G 2 (0 



C_ 

„4 ! 



(18) 



where C is a constant of integration. 

Variation with respect to the lapse function N and 
setting N = 1 at the end, leads to the first Fricdmann 
equation 



3M Pl H 2 



+ p r - 3m 2 M Pl Gi 



1 
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ZG', 



(19) 



where primes in G"s denote derivatives with respect to 
their argument £. Similarly, variation with respect to a 
provides the dilaton evolution equation 

(a 3 a) + UI Pl m 2 [Fi(0 + afF 2 (d = 0, (20) 



a 3 Ndt 
where 

Fi(0 = 3(1 + ja 3 + a 4 )£ - 2(1 + ^a 3 + 3a 4 )£ 2 

+ ^(a 3 +4a 4 K 3 (21) 

F 2 (0 = (1 + ^a 3 + Q- 4 )C - 2(1 + ^a 3 + 3a 4 )£ 2 



-(a 3 +4a 4 )£ 3 -4a 4 £ 4 



(22) 



Finally, the evolution equations close by consider- 
ing the continuity equations for matter and radiation, 
namely p m + 3H(p m +p m ) = and p r + 3H(p r +p r ) = 
respectively. 

Before proceeding to the detailed analysis of the model, 
let us make a few comments. Firstly, as it was first no- 
ticed for a similar model in [671 ] , the Fricdmann equation 
(|19[) can exhibit singularities when the denominator goes 
to zero. However, as we will show in the next section, 
these singularities separate the phase space in discon- 
nected branches, that is the universe cannot evolve from 
one to the other. 

As a second remark we mention that in the case where 
C^O, in the constraint equation (|18[) it appears a func- 
tion of 1/a 4 which plays, in the early Universe, a role 
similar to "Dark Matter" or "Dark Radiation". In par- 
ticular, when a 4 ^ 0, in the early Universe (a — > 0) rela- 
tion (|r7|) gives £ — > oo. This implies G 2 — — «4£ 4 = C/a 4 
leading the Friedmann equation (fj~9"|) to be 



3M Pl H 2 



3m 2 M 2 



pi 



. « 3 



C\3/4 



OS4, 



that is an effective dark-matter term appears. Similarly, 
when a 4 = we acquire 

3 (l - ^jM Pl H 2 ~ p m + p r + m 2 M^Ca- 4 , (24) 

that is we get an effective dark radiation term. 

Finally, we mention that in the particular case C = 0, 
(|18p leads to a constant £, and the Friedmann equation 
([H?]) becomes 



3A/|,,iJ 2 = p m + p r -a A + H 2 



(25) 



where a\, j3 are constants. It is interesting to notice that 
this relation appears as a particular case of Holographic 
models where the IR cut-off is fixed at the Hubble scale. 
In fact in the aforementioned theories the number of de- 
grees of freedom in a finite volume should be finite and 
related to the area of its boundary. This gives an upper 
bound on the entropy contained in the visible Universe. 
Following an idea suggested by [82], [83( where a long dis- 
tance cut-off is related to a short distance cut-off, in J84J 
it was suggested to take the largest distance as the Hub- 
ble scale or the event horizon. Later these models were 
dubbed holographic dark energy [85I J . 

All these particular cases reveal the richness and the 
capabilities of the scenario. In what follows we shall cast 
the evolution equations in the autonomous form and in- 
vestigate their implications for late time cosmology. 



III. DYNAMICAL BEHAVIOR 

In this section we perform a dynamical analysis of the 
scenario at hand, which allows us to bypass the com- 
plexities of the equations and obtain information for the 



(23) 



late-time, asymptotic, cosmological behavior [86M88J . We 
mention here that in [78| the authors investigated the cos- 
mological implications in a specific setting allowing them 
to obtain simple analytical solutions. However, in this 
work we are interested in the full dynamical investiga- 
tions of evolution equations. 



A. Phase-space analysis 

In order to perform a phase-space analysis, we first 
transform the involved cosmological equations into their 
autonomous form introducing suitable auxiliary vari- 
ables. We extract the critical points of the dynamical 
system under consideration, perturb the system around 
them and analyze their stability by examining the eigen- 
values of the corresponding perturbation matrix. 



Let us firs introduce the density parameters 

Pni 



^ t"m 



3M Pl H 2 

Pr 

3M' Pl H 2 



O 2^1 



(26) 
(27) 

(28) 
(29) 



with which we can re-write the Friedmann equation (|19p 
as 



- -i, 



il r -\- ill 



0, = l. 



(30) 



Thus, we can now transform the above cosmological sys- 
tem into its autonomous form, using only the dimension- 
less variables ft r , $1a and £, while (|30|) will be used to 
eliminate f2 m . Doing so we obtain 



d£l r 
din a 
d$l A 
din a 

din a 



H 2+ w) 




a A (s 
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H 
H 2 


G 2 






G' 2 ' 







(31) 
(32) 
(33) 



where the combination -Mj can be acquired differentiating 
the first Friedmann equation (fT9f as 



H 

H 2 



%J A L fii J_ £j _; - j 



12-22. 



Gi 



5r^A + f^(i-4^ 



CGJ, 



(34) 



A comment about the role of dilaton in the scenario un- 
der consideration is in order. The dilaton role is twofold, 
first, it has an effect through its impact on the graviton 
mass, since in the above equations m appears multiplied 
by a function of £, that is of the exponential of the dila- 
ton field. This is what we call fl\. Secondly, the dilaton 
has an effect through its kinetic term in (fT4| , which is of 
course proportional to the parameter to. This is what we 
call fl a . Therefore, strictly speaking, in the present sce- 
nario the effective dark energy sector includes both $7a 
and f^, since these two contributions cannot be distin- 
guished observationally (we write them separately only 
for clarity). Thus: 



VLde = Qa + O ct = 1 — O r 






(35) 



Finally, note that in terms of the auxiliary variables 
f2 r , SI a and £, we can express the physically interest- 
ing observables such as the total (effective) equation-of- 
state parameter w e ff, the deceleration parameter q and 
the dark energy equation-of-statc parameter wde respec- 
tively as 



W eff 



-1 



2 H 
3H 2 



(36) 



q = -l- 



H 

H 2 



WDE 



W eff 



1 3 

= O + oWeff- 



UU'P u uf 



sia +n„ 



(37) 



(38) 



where the combination H / H 2 is given by (jM)) . Lastly, in 
the following for simplicity we focus on the usual cases 
of standard dust matter w m = and standard radiation 
w r = 1/3. 

The critical points of the above autonomous system are 
extracted setting the left hand sides of equations (f3"T|) - 
(|33p to zero. In particular, equation (|33[) implies that 
at the critical points G2(£) = 0, which according to (fT6|) 
gives £ = 0, £ = 1 and 



£t = i 



3^3 ± yjQa 2 — 64a 4 

8c*4 



(39) 



In order for the critical points to be physical they have to 
possess < fi r < 1, < SIa+Oo- < 1, < £ and of course 
£ £ 1Z. Now, in order for £± £7?we require 9a§ — 64a4 > 
while < £ leads to the necessary conditions for a's 
using ([55]). 

The real and physically meaningful critical points 
(fi r ,fiA,£) of the autonomous system (j3"Tj) -([53" |) are pre- 
sented in Table U along with their existence conditions. 
We mention here that there is an additional existence 
condition, namely the obtained H 2 from (fT9"|) must be 
always positive, a condition which constrains the allowed 
as and 0:4 values. We will come back to this later on. 

For each critical point we calculate the 3x3 perturba- 
tion matrix of the corresponding linearized perturbation 
equations, and examining the sign of the real part of its 
eigenvalues we determine the type and stability of this 
point (positive real parts correspond to unstable point, 
negative real parts correspond to stable point, and real 
parts with different signs correspond to saddle point). 
The eigenvalues and the stability results are summarized 
in Table |U 

Finally, for completeness, in the same Table we present 
the corresponding values of the rest density parameters 
Cl m and Q a , as well as the calculated values for the ob- 
servables such as the deceleration parameter q, the total 
cquation-of-state parameter w e g, and the dark energy 
equation-of-statc parameter Wde- These results are in 
agreement with the general analysis of nonlinear massive 
gravity of [831 ■ 

Before proceeding to the discussion of the physical im- 
plications of these results, we make a remark concerning 
the structure of the phase space. As we mentioned in the 
Introduction, the Friedmann equation (|19[) exhibits sin- 
gularities when the denominator goes to zero. In particu- 
lar, in the scenario at hand we obtain various singularities 

localized at the roots of 1 — ^ ( 1—4™- ) , which, for every 

parameter w, is a function of the rescaled field- variable £ 
only. Since ui must be relatively small in order to satisfy 
observations (this will be confirmed later on), assuming 
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1 - 3w/2 
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< w < 2/3 


Saddle point 


-4, 3, -1 


B 





1 - 3w/2 








3w/2 
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-1 


-1 


< w < 2/3 


Attractor 


-4, -4, -3 
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1 - 3w/2 











3w/2 
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Saddle point 


-4, 4, 1 
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l-w/6 
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< w < 6 


Attractor 


-4, -1, -1 


E± 
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-4,-1, 3 
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1/3 
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-4,4,1 
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0< w <6, <£±, £± G^ 
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-4,-4,-3 
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n r 


1 — w/6 — fi r 
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w/6 
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1/3 


1/3 


< w < 6 


Saddle 


-4,1,0 



TABLE I: The real and physically meaningful critical points of the autonomous system (|31|) -(|33 |) . their existence conditions, 
the eigenvalues of the perturbation matrix, and the deduced stability conditions. We also present the corresponding values of 
the various density parameters, and the values of the observables: deceleration parameter q, total equation-of-state parameter 
w e ff and dark energy equation-of-state parameter wde- 



w< 1 the zeros of the previous function nearly coincide 
with the roots of G' 2 ((,)- But as we mentioned above, the 
critical points correspond to G%(£) = 0. Thus, using the 
Rolle's theorem, we conclude that each cosmological late- 
time solution (G2(£) = 0) and the corresponding region, 
is disconnected from the others by a singularity (G' 2 = 0). 
In other words, during its evolution, the universe cannot 
transit from one of these regions to another. 



namical systems [90, |9l| , with each critical point attract- 
ing orbits initially starting inside its basin of attraction, 
in the present scenario the case is simpler, since these 
points usually lie in the disconnected regions we men- 
tioned in the previous subsection, separated by the £ val- 
ues that make the denominator of (fT^) zero, and thus to 
Big-Rip-type divergences [92H95I ]. Finally, note that for 
uj > 6 there are not stable late-time solutions, that is the 
corresponding cosmology is not interesting. 



B. Physical implications 

In the previous subsection we performed a complete 
phase-space analysis of the scenario at hand, we ex- 
tracted the late-time stable solutions (attractors) and we 
calculated the corresponding observables. Here we dis- 
cuss their cosmological implications. 

As we observe from Table HI there exist four stable crit- 
ical points, that can attract the universe at late times, 
namely B, D and the two points G±. Points B and 
G± correspond to a dark-energy dominated, accelerat- 
ing universe, where dark energy behaves like cosmologi- 
cal constant (de Sitter solutions). These features make 
these points very good candidates for the description of 
late-time universe. On the other hand, point D corre- 
sponds to a non-accelerating universe, with dark energy 
and matter density parameters of the same order (that 
is it can alleviate the coincidence problem) , but the dark 
energy has a stiff equation of state. Therefore, this point 
is disfavored by observations. 

Apart from the above stable late-time solutions, the 
scenario of quasi-Dilaton non-linear massive gravity 
possesses the non-accelerating, stiff dark-energy points 
A,C,E±,F±, and the critical curve /. These points are 
saddle (/ is non-hyperbolic curve with saddle behavior) 
and thus they cannot be the late time states of the uni- 
verse, and moreover their cosmological features are not 
favored by observations. 

Let us make a comment here on the appearance of 
more than one stable critical points, which for some pa- 
rameter choices can exist simultaneously. Although the 
bi-stability and multi-stability is usual in complicated dy- 



-de Sitter 
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^t-n, 



FIG. 1: Trajectories in the ilDE-Qm plane of the cosmological 
scenario of quasi-Dilaton non-linear massive gravity, for the 
parameter choice as = \,o.a = 1 and w = 0.01. In this spe- 
cific example the universe at late times is led to the de Sitter 
solution. 



Finally, as we previously mentioned, the 013-04 param- 
eter subspacc is constrained by the condition the ob- 
tained H 2 from (fl"9|) to be always positive. The cor- 
responding constraints are complicated and we do not 
present them here, however they do get a very simple 
form when the universe evolves towards its interesting 
stable late-time states B and G± (where f2 r and fi m be- 
come smaller and smaller). In particular, for the point B 



we obtain 2 + a.3 + a.4 < 0, for the point G+ we acquire 

2 

«3 > 0, < ay < -# and for the point G_ we obtain 

2 

«3 < and < a4 < ^-. 

In order to present the cosmological behavior more 
transparently, we numerically evolve the autonomous sys- 
tem pTj) - ([55| for the the parameter choice 0:3 = 1,04 = 1 
and w = 0.01, and in Fig. Q]we depict the corresponding 
phase-space behavior in the floE-^m plane. 



IV. COSMOLOGICAL EVOLUTION 

In the previous section we performed a dynamical anal- 
ysis of quasi-Dilaton non-linear massive gravity, that is 
we focused on its asymptotic behavior at late times, that 
is on solutions that are going to attract the universe in- 
dependently of the initial conditions and of the specific 
cosmological evolution towards them. In this section we 
investigate the behavior of the universe at intermediate 
times, which obviously does depend on the initial condi- 
tions, but it can be very interesting. In particular, we 
are interested in obtaining a evolution of the universe in 
agreement with the observed epoch sequence. Further- 
more, we examine the possibility of obtaining bouncing 
behavior. In the following subsections we discuss these 
two cosmological behaviors separately. 



A. Epoch sequence and late-time acceleration 

Let us first examine the universe evolution, focusing 
on the various density parameters and the dark-energy 
equation of state. Observing the Fricdmann equation 
([12]), as well as the evolution equations for the density 
parameters and for £ ([3Tj) - (|33|) . we deduce that the ob- 
served post-inflationary thermal history of the universe 
can be easily obtained by suitably choosing the model 
parameters. 

In order to present this behavior more transparently, 
we numerically evolve the system for 0J3 = 1, 0:4 = 0.115 
and ui = 0.01, starting with initial conditions correspond- 
ing to Q r f=a 1, and imposing the current values to be 
O m « 0.28 and fl DE ss 0.72. In Fig. [2] we depict the 
resulting evolution for the density parameters using, in- 
stead of the scale factor, the redshift z as the indepen- 
dent variable (1 + z = a^/a with ao = 1 the present 
scale- factor value). As we observe, the thermal history 
of the universe can be reproduced, namely we obtain the 
successive sequence of radiation, matter and dark energy 
epochs. 

Similarly, in Fig. [3] we present the corresponding evo- 
lution of the total equation-of-state parameter w e ff as 
well as of the dark-energy one wde- As we observe, 
dark energy starts from a dust-like behavior, resulting 
in a cosmological-constant-likc one at the current uni- 
verse, driving the observed acceleration. Note that, as 
expected, in the future, where dark energy completely 



1.0 


;~a 


1 » *H1 >v 


s — ' 


0.8 


; \ 




v / flr 


0.6 


- 




\ / 


0.4 


- 




/ \ 


0.2 
0.0 


J 


\ ^DE 


' \ 








-1 



2 3 4 

Log 10 [1 + z] 



FIG. 2: Cosmological evolution of the radiation (Cl r ), matter 
(Vim) and dark energy (T2.de = ^2a + &■&) density parameters, 
as a function of the redshift z, for the parameter choice 03 = 
1, 04 = 0.115 and ui = 0.01. In this specific example the 
observed universe epoch sequence is reproduced. 
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FIG. 3: Cosmological evolution of the total (w e ff) and of the 
dark-energy (wde) equation-of-state parameters, as a func- 
tion of the redshift z, for the parameter choice az — 1, a& = 
0.115 and uj = 0.01. In this specific example, dark energy 
starts from a dust-like behavior, resulting in a cosmological- 
constant-like one at the current epoch, which drives the ob- 
served universe acceleration. 



dominates, both w e ff and wde coincide at the cosmo- 
logical constant value. 



B. Bouncing and turnaround solutions 

In this subsection we examine the bounce realization 
in the scenario at hand. Bouncing cosmologies in gen- 
eral have gained significant interest, since they offer an 
alternative paradigm free of the "Big-Bang singularity" , 
as well as of the horizon, flatness and monopolc problems 
[9a |. The basic condition for the bounce is the Hubble 
parameter to change sign from negative to positive, while 



its time-derivative to be positive. Additionally, when 
the Hubble parameter changes from positive to nega- 
tive, while its time-derivative is negative, we have the 
realization of a cosmological turnaround. Such a behav- 
ior implies the Null Energy Condition (NEC) violation 
around the bounce (or the turnaround), which is non- 
trivial in the context of General Relativity I97J| . leading 
in general to ghost degrees of freedom [98|, |99[. How- 
ever, a safe bounce evolution, free of ghosts and instabil- 
ities, is po ssible in various modified gravity constructions 
jlOOl Il0l| . Clearly, observing the Friedmann equation 
(HHJ) of the present scenario of quasi-Dilaton non-linear 
massive gravity, we deduce that the basic bounce (or 
turnaround) realization condition can be easily fulfilled. 
In order to quantify the above discussion, and assum- 
ing standard dust matter and standard radiation, we 
rewrite the Friedmann equation (|19[) in the following 
form: 



Hi 



a, 
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(40) 



where £l m o and £l r ,o are the values of the correspond- 
ing density parameters at present (ao = 1), and Hq the 
current Hubble parameter. Furthermore, since the con- 
straint equation (fT5|) provides £ as a function of the scale 
factor, we can finally write (|40|) as 



IP 
HI 



V(a), 



(41) 



with C, (33, (34, u, m and £l m o,Sl r .o as parameters. In 
summary, examining the behavior of the above known 
V{a) 1 we can find possible bouncing points at a = a B 
(when V(a B ) = and H^ 1 j- t ^/V{a B ) > 0), or possi- 
ble turnaround points at a = ax (when V(o,t) = and 
H ~ 1 j^y / V(aT) < 0). This procedure is straightforward, 
and thus we deduce that both bounce and turnaround are 
possible in the scenario of quasi-Dilaton non-linear mas- 
sive gravity. Although one can apply it in full generality, 
obtaining exact results, since the involved expressions for 
a B and ax are lengthy in the following we explicitly apply 
it in a specific simple example, where simple expressions 
can be easily obtained. 

First of all, since a bounce occurs in the early universe, 
that is at small scale factors (a <C ao = 1), from (fT5|) we 
deduce that at the bounce £ S> 1. On the other hand, 
since a turnaround occurs in the late universe, that is at 
o > ao - 1) from (|18p wc deduce that at the turnaround 
G 2 (0 = 0. 

Let us for simplicity investigate the simplest case a% = 
«4 = 0. At the bounce region from (fT5| we obtain G\ « 
i 2 and G 2 « -<f = C/a 4 , which implies that C < 0. 
Therefore, from (|4"0"|) we deduce that 
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with 



and therefore 
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From the above relations we easily acquire that H(a B ) = 
and moreover for cj < 6, which is always the case in in- 
teresting cosmology as we discussed in section IIII1 we 
have H(clb) > 0. Thus, we immediately see that the 
bounce is obtained at a — a B . Finally, note that ob- 
viously we require a B > 0, which gives an additional 
constraint on the parameters, namely Cm 2 /H^ < — il r .Q. 
Let us now see whether a turnaround can occur in this 
simplest case «3 = a 4 = 0. For a >• 1, that is for G2 = 0, 
we have £ = or £ = 1. Obviously, since £ varies, the 
above conditions are quite difficult to be obtained, and 
thus contrary to the bounce, the turnaround is harder 
to be acquired. In the £ = case from (|15[1 we acquire 
G\ = 2, and therefore from pi)]) we obtain 
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where for simplicity we neglected the radiation term since 
at large scale factors the matter term will dominate in 
general. Thus, 
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From these relations we deduce that H(ax) = 0, while 
H(o.t) < for u; < 2/3, which is the turnaround condi- 
tion (note that in this case or > always). Concerning 
the second case £ = 1, from (fl~5j) we acquire G\ =0, 
and thus H 2 oc fl m oa~ 3 + £7 r! oa -4 which has no roots, 
and hence we deduce that in this case the turnaround is 
impossible. 

Finally, note that for w < 2/3 and Cm 2 / 'Hq < —il r fi, 
in the above simplest example we obtain both a bounce 
and a turnaround, which is the realization of cyclic cos- 
mology [9g|. In Fig. 21 wc present such a behavior, by 
numerically evolving the exact cosmological equations for 
parameter choices that satisfy the above conditions. 

One can straightforwardly perform the above analysis 
for general a.3 and 04, extracting the corresponding con- 
ditions for a bounce and turnaround, or for their simul- 
taneous realization, that is for cyclic cosmology. In sum- 
mary, as we can see, in the scenario of quasi-Dilaton non- 
linear massive gravity the above alternative cosmological 
evolutions can be easily obtained, in agreement with the 
general case of extended non-linear massive gravity [77| . 
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FIG. 4: The rescaled Hubble function h = H/Hq as a function 
of a/ao, for Q3 = 04 = 0, fi m = 0.3, fi r = 1CP 4 , C — 
—1.1 x 10~ 4 and m = Hq. In this case we obtain two evolution 
branches, one corresponding to a cyclic Universe (blue line) 
and a second branch corresponding to a Universe with a past 
bounce (red dashed curve). The dotted line corresponds to 
A.CDM scenario with the same cosmological parameters. 



V. OBSERVATIONAL CONSTRAINTS 

Having analyzed the cosmological behavior in the sce- 
nario of quasi-Dilaton non-linear massive gravity, we 
now proceed to investigate the observational constraints 
on the model parameters. We use observational data 
from Type la Supernovae (SNIa), Baryon Acoustic Os- 
cillations (BAO), and Cosmic Microwave Background 
(CMB), and as usual we work in the standard units 
suitable for observational comparisons, namely setting 
8ttG = c= 1. 

In order to perform the analysis, and assuming stan- 
dard dust matter and standard radiation (w r = 1/3), 
we re-write the Friedmann equation (|19[) in the form 
of (|4TJ)l . and we use the redshift z instead of the scale 
factor. Thus, the scenario at hand has the parameters 
C, 03, 04, u, m and f2„ ,,n or ^deo ( we fix H by its 7-year 
WMAP best-fit value |103|). The details of the combined 
SNIa+CMB+BAO analysis are given in the Appendix, 
and here wc present the constructed likelihood contours, 
discussing their structure. 

Let us first examine the constraints on 03 and 04. In 
Fig. [5]wc depict the la and 2a 0:3 — Q m o and 04 — Q m o 
contours respectively. As we observe, the parameters 03 
and 04 are not constrained by the observational data, 
that is they only slightly affect the cosmological evolu- 
tion. This behavior verifies the results obtained in section 
IIII1 where the stable late-time solutions were found to be 
independent of 013,04, as well as those obtained in section 
IIVI where at intermediate times different values of 03,04 
lead to nearly the same cosmology. This is in agreement 
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FIG. 5: Top: The la (dark) and 2a (light) likelihood contours 
in the 03 — Q. m o plane, using SNIa, BAO and CMB data, for 
Ola = 1. Bottom: The la (dark) and 2a (light) likelihood 
contours in the an — Q. m o plane using SNIa, BAO and CMB 
data, for 03 — 1. Both graphs are for ui — 0.01. 



with the analysis of |57j for the simple non-linear mas- 
sive gravity. Finally, concerning $7 m o we obtain a best fit 
value approximately at 0.284 for both 03 and 04. 

In order to extract the constraints on the quasi-dilaton 
parameter uj, and knowing the above results on the in- 
significant role of 03 and 04, we fix them to 03 = 1,04 = 
0.115 (alternatively we could consider them as free pa- 
rameters, marginalizing with fiat priors, but the results 
would be almost the same). In Fig. [S]we depict the la 
and 2cr likelihood contours in the uj — Q m Q plane. As 
we observe u> is significantly constrained by observations, 
with its best-fit value being ui = 0.0878. Note also that 
negative uj values are also allowed, leading the dilaton 
to behave as a phantom field. This strong constraining 
behavior verifies the results of the dynamical analysis of 
section HIT1 where in the case of interesting late-time cos- 
mology ui was bounded. Finally, we mention that apart 
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FIG. 6: The la (dark) and 2a (light) likelihood contours in the 
u — QmO plane using SNIa, BAO and CMB data, for 03 = 1 
and «4 = 0.115. 



from the above late-time observational data, we could 



use Big Bang Nucleosynthesis arguments |103| in order 
to impose constraints on £Ide, that is to w, however these 
would be weaker than the ones obtained above. 

Finally, we examine the constraints on the graviton 
mass 77i, and in Fig. [7] we depict the corresponding like- 
lihood contours in relation to fl m o- From this figure we 
observe that the ratio tu/Hq is highly constrained, and 
in particular we verify the expected result m ~ Hq, al- 
though with m being a bit smaller than H . In par- 
ticular, the best-fit values are iti/Hq = 0.44775 and 
il m0 = 0.2794. 
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FIG. 7: The la (dark) and 2a (light) likelihood contours in 
the m/Ho — fi m o plane using SNIa, BAO and CMB data, for 
«3 = 1 and a.4, = 0.115. 



VI. CONCLUSIONS 

In this work wc investigated the cosmological dynamics 
of the quasi-Dilaton non-linear massive gravity. In such a 
gravitational construction, one uses the usual non-linear 
massive gravity [y, |7j , which is free of ghosts and couples 
it with a dilaton field [7g. Apart from the matter and 
radiation sectors in this scenario we obtain an effective 
dark energy, constituted from the massive graviton and 
the dilaton contributions. The corresponding cosmologi- 
cal behaviors proves to be very interesting. 

First we performed a detailed dynamical analysis, and 
we extracted the stable late-time solutions, that is solu- 
tions that always attract the universe at late times. Fur- 
thermore, in each of these states we have calculated the 
corresponding obscrvables, such as the dark-energy and 
matter density parameters, the dark-energy equation-of- 
state parameter, and the deceleration parameter. We 
showed that the scenario at hand gives rise to a dark- 
energy dominated accelerating universe where dark en- 
ergy behaves like cosmological constant. 

Leaving aside the good asymptotic behavior at late 
times, the cosmological behavior of quasi-Dilaton non- 
linear massive gravity has also interesting intermediate 
epochs a la thermal history. We have shown that the 
observed history of the universe can be easily obtained 
in a parameter region, that is the expected sequence of 
radiation, matter, and dark energy epochs. As a bonus, 
we found conditions on the parameters that might lead 
to bouncing, turnaround or cyclic behavior which is an 
interesting feature of the scenario under consideration. 

Last but not least, we have imposed observational con- 
straints on the model parameters. We carried out model- 
fitting using data from Type la Supernovae (SNIa), 
Baryon Acoustic Oscillations (BAO), and Cosmic Mi- 
crowave Background (CMB) observations. Construct- 
ing the corresponding contour plots we deduced that the 
usual 03 and a^ parameters of non-linear massive gravity 
arc not constrained by current data at the background 
level. However, the dilaton-coupling parameter is sig- 
nificantly constrained at a narrow window around zero. 
Finally, as expected, the graviton mass is found to be of 
the same order of the present value of the Hubble param- 
eter. 

From the aforesaid, we find that cosmology of the 
quasi-Dilaton non-linear massive gravity is interesting 
and is in agreement with the observed universe behav- 
ior. In our opinion, it can be a good candidate for the 
description of universe. However, in order to obtain ad- 
ditional indications of the validity and consistency with 
observations, of the scenario at hand, apart from the 
background analysis that was the base of this work, we 
should proceed to a detailed examination of the pertur- 
bations. In particular, one would need to investigate the 
perturbation-related observables such as the spectral in- 
dex, the tcnsor-to-scalar ratio, and the non-gaussianity 
estimators. Such an analysis lies beyond the goal of the 
present work and it is left for future investigations. 
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Appendix: Observational data and constraints 

Here we briefly review the sources of observational con- 
straints used in this manuscript, namely Type la Super- 
novae constraints, Baryon Acoustic Oscillations (BAO) 
and Cosmic Microwave Background (CMB). The total 
X 2 defined as 



X 



XsN 



Xbao 



XCMB j 



(A.l) 



where the individual xf f° r every data set is calculated 
as follows. 

a. Type la Supernovae constraints 

In order to incorporate T ype la constraints we use the 
Union2.1 data compilation 104j of 580 data points. The 
relevant observable is the distance modulus /i which is de- 
fined as \i — m — M = 5 log Dl + Ho, where m and M are 
the apparent and absolute magnitudes of the Supernovae, 
Dl{z) is the luminosity distance Dl{z) = (1+z) J Q Z jjtt\ 

and no = 5 log I -jf— I + 25 is a nuisance parameter that 
should be marginalized. The corresponding x 2 writes as 

xsn{ij>o,v) = 2^ _. ( -\2 ' ( A - 2 ) 

i=1 °tl\z l ) 

where Hobs denotes the observed distance modulus while 
fith the theoretical one, and a^ is the uncertainty in the 



distance modulus. Additionally, 9 denotes any parameter 
of the spe cific model at hand. Finally, marginalizing /j,q 
following [l05l | we obtain 



X 2 SN {e) = A{6) - 



B(6) 2 



(A.3) 



with 



^ o~^f 

IMh(zi,Ho = 0,0) - Hobs{zi) 



8=1 

580 



b(9) = y: 



580 



0>(z») 2 



C{0)=Y— —. 

^ °y.(*) a 

b. Baryon Acoustic Oscillation constraints 
We use BAO data from [l06l4T09J . that is of ■ 



(A.4) 
(A.5) 
(A.6) 



<Ia(z*) 



D v (Z BAO )' 

where z* is the decoupling time given by z+ w 1091, the 
co-moving angular-diameter distance oIa{z) = J Q jttt) 

and Dy(z) = (dA(z) 2 jjf^) is the dilation scale |llOJ |. 
th required data are depicted in Table 1771 

In order to calculate Xmao f° r BAO data we follow the 
procedure described in (lllj . where it is defined as, 
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BAO^BAO 



)XbAO, 
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with 
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(Ja(z*) 
£>y(0.106) 
d A (z*) _ 
Dvf0.2) 

dA(Zt) 

Dv(0.35) 
dAJz*) 

D V (0A4) 
dA(z*) 
Dv(0.6) 
(Ja(z*) 

D v (0.73) 



- 30.95 \ 

- 17.55 

- 10.11 
-8.44 

6.69 
-5.45 



(A- 



and the inverse covariance matrix reads as 



J 



c- 



0.48435 


-0.101383 


-0.164945 


-0.0305703 


-0.097874 


-0.106738 \ 


-0.101383 


3.2882 


-2.45497 


-0.0787898 


-0.252254 


-0.2751 


-0.164945 


-2.45499 


9.55916 


-0.128187 


-0.410404 


-0.447574 


-0.0305703 


-0.0787898 


-0.128187 


2.78728 


-2.75632 


1.16437 


-0.097874 


-0.252254 


-0.410404 


-2.75632 


14.9245 


-7.32441 


-0.106738 


-0.2751 


-0.447574 


1.16437 


-7.32441 


14.5022 / 



(A.9) 



c. CMB constraints 
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ZBAO 


0.106 0.2 0.35 0.44 0.6 0.73 


dAM 

D \( Z BA0) 


30.95 ± 1.46 


17.55 ±0.60 


10.11 ±0.37 


8.44 ± 0.67 


6.69 ± 0.33 


5.45 ±0.31 



TABLE II: Values of 



D v( Z BAo) 



for different values of zbao- 



We use the CMB shift parameter 

/"1089 dz i 

R = H \/ f2 m0 / 
Jo 



(A.10) 



/o H(z') 

following J102J . The corresponding Xcmb ^ s defined as, 

(R(0) - i? ) 2 



Xcmb(0) 



(A.ll) 



with R = 1.725 ± 0.018 and R(6) [10 
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